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Abstract 

We study the interplay between flavor symmetries and leptogenesis in the case 
when the scale of flavor symmetry breaking is higher than the scale at which lepton 
number is violated. We show that when the heavy Majorana neutrinos belong to an 
irreducible representation of the flavor group, all the leptogenesis CP asymmetries 
vanish in the limit of exact symmetry. In the case of reducible representations we 
identify a general condition that, if satisfied, guarantees the same result. We then fo- 
cus on the case of a model in which an flavor symmetry yields a drastic reduction 
in the number of free parameters, implying that at leading order several quantities 
are only a function of the lightest neutrino mass mi, which in turn is strongly con- 
strained. For normal ordering (NO) we find mi ~ (0.0044 0.0056) eV while for 
inverted ordering (10) mi > 0.017 eV. For the Oi^2/3 decay parameter this yields 
|mee| ~ (0.006 ^ 0.007) eV (NO) and \mee\ > 0.017eV (10). We show that the lep- 
togenesis CP asymmetries only depend on m^, on a single non- hierarchical Yukawa 
coupling y, and on two parameters that quantify the flavor symmetry breaking ef- 
fects, and we argue that the unflavored regime for leptogenesis is strongly preferred 
in our model, thus realizing a rather predictive scenario. Performing a calculation 
of the matter-antimatter asymmetry we find that for NO the observed value is eas- 
ily reproduced for natural values of the symmetry breaking parameters. For 10 
successful leptogenesis is possible for a limited choice of the parameters implying 
rather large reheating temperatures Tj-eh > 5 X 1013 GeV. 



1 Introduction 



Perhaps the simplest explanation of the smallness of neutrino masses is in terms of the 
largeness of the scale at which B — L is violated, through the well-known seesaw mech- 
anism [1] that, in the three neutrino framework, also implies CP violation in the lepton 
sector. Barring very special cancelations such a large scale is expected to be well above 
the weak scale, possibly near or just below the grand unified scale. It is an intriguing 
coincidence that also baryogenesis, if occurring at temperatures well above the weak scale, 
needs B — L violation, since any initial B + L asymmetry would be erased by sphaleron 
interactions in the subsequent evolution of the universe. The seesaw mechanism and 
baryogenesis elegantly combine in leptogenesis [2] where the observed baryon asymmetry 
is produced by the out-of equilibrium, CP and B — L violating decays of the right-handed 
(RH) neutrinos. It is also quite remarkable that, at least in its simplest implementation, 
leptogenesis requires light neutrino masses below the eV scale [3l HI [5] , in a range which is 
fully compatible with other experimental constraints. Unfortunately, it is difficult to pro- 
mote this elegant picture into a testable theory, due to the large number of independent 
parameters of any seesaw model. 

For this reason it is of great interest to consider more constrained frameworks, where, 
by restricting the number of relevant parameters, the observed baryon asymmetry can 
be related to other observable quantities. Examples of such frameworks are models of 
lepton masses where the total number of parameters is reduced by the presence of a flavor 
symmetry. There is a vast variety of models of this type: the flavor symmetry can be 
global or local, continuous or discrete and there are many viable candidate groups. A 
common feature of these models is that the symmetry is always broken, to provide a 
realistic description of the lepton masses and mixing angles. The breaking can be either 
explicit or spontaneous and is generally controlled by a set of dimensionless symmetry 
breaking quantities rj. In general any such model can be thought of as an expansion 
in powers of t]. There is a leading order approximation, corresponding to the limit of 
exact flavor symmetry, plus small corrections proportional to rj and its positive powers. 
If T) are small, it makes sense to neglect high powers of rj and the model predictions will 
effectively depend on relatively few parameters. Many models of lepton masses based on 
flavor symmetries reproduce the most relevant features of the lepton spectrum. The main 
idea is to describe small dimensionless observable quantities, such as charged lepton mass 
ratios, ^13, 623— it /A, in terms of some positive power of the symmetry breaking parameters 
rj. In the context of leptogenesis, given the extreme smallness of the baryon asymmetry, 
it is quite natural to expect that the CP asymmetries in the RH neutrino decays are also 
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suppressed by some power of 77, opening the interesting possibility of relating, through 
specific flavor symmetries and flavor symmetry breaking patterns, the baryon asymmetry 
to others low-energy observables. Indeed if these CP asymmetries were of order 77°, we 
would need an abnormal suppression coming either from a dilution factor from photon 
production or from the wash-out factor. Typically we expect a dilution factor of order 
10~^ and a wash-out factor in the range 10^^ ^ 10~^, which favor CP asymmetries around 
10~6 ^ 10"^ 

We think that establishing the dependence of the CP asymmetries on rj could be of 
great interest for model building and for phenomenology. In this paper we collect a number 
of observations that can be useful to this scope. In the first part of this work, section 2, 
we consider a completely general framework and we derive conditions for the vanishing 
of the total CP asymmetries based on general group theoretical arguments. The main 
result of this part is the proof that when the heavy Majorana neutrinos are assigned to an 
irreducible representation of the flavor group, all the leptogenesis CP asymmetries vanish 
in the limit of exact symmetry. A general condition for the vanishing of the leptogenesis 
CP asymmetries that hold for different types of flavor symmetries is also given. 

In the second part of the paper we work out in detail the baryon asymmetry predicted 
by a model of lepton masses based on the A4 discrete symmetry group [6], originally 
built to reproduce the tri-bimaximal lepton mixing (TBM), a good approximation of the 
existing data [7]. In such a model there are strong indications that leptogenesis can be 
successfully realized, since the leptogenesis CP asymmetries are of order t]'^ [8] and the t] 
parameters are of order 10~^. However, an estimate of the washout effects is still lacking 
and, given the small number of independent parameters, it is not a priori guaranteed that 
the baryon asymmetry can be reproduced. 

In section 3 we review the main features of the A4 model. We show how this model 
is surprisingly predictive: if the neutrino ordering is normal, the neutrino mass spectrum 
and the Oz/2/? decay parameter | determined rather precisely. If the ordering is 

inverted, interesting lower limits on both quantities are obtained. In section 4 we confront 
the A4 model with leptogenesis. We first compute the leptogenesis CP asymmetries at 
sub- leading order. Next, both for the normal and inverted ordering cases, we estimate the 
value of the cosmic baryon asymmetry resulting from leptogenesis. It is quite surprising 
that in the case of normal hierarchy the right values for the CP asymmetries are obtained 
despite the tight constraints found on the involved parameters, namely on the phases and 
on the absolute neutrino mass scale. In the case of inverted ordering successful leptogenesis 
is possible but it relies on a compensation between very large CP asymmetries \ei\ ~ 10~^ 
and strong wash-out from AL = 2 processes. This compensation occurs for very large 
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RH neutrino masses ~ 10^ GeV implying in turn values of the reheating temperature not 
lower than ~ 5x 10^^ GeV for quasi-degenerate light neutrinos, only marginally compatible 
with cosmological observations. Finally, in section 5 we collect our main results and draw 
the conclusions. 



2 Flavor symmetries and leptogenesis CP asymme- 
tries 

In this section we discuss the constraints on the leptogenesis CP asymmetries implied by 
different classes of flavor symmetries. We assume that the flavor symmetry is broken at a 
scale above the leptogenesis scale, that is, leptogenesis occurs in the standard framework of 
the SM plus three heavy Majorana seesaw neutrinos [1], without involving in its dynamics 
any additional stateQ The most interesting situations occur when the leptogenesis CP 
asymmetries vanish in the limit of exact flavor symmetry. Since in several models the 
breaking of the symmetry is related in a clear way to parameters that are measurable at 
low energy (like e.g. 6*13 or 7r/4 — 6^23) our results open up the possibility of constraining 
- within specific flavor models - the leptogenesis CP asymmetries (that are the crucial 
high energy parameters) from low energy experiments. For definiteness we refer to the 
non-supersymmetric case but all our results hold in the supersymmetric case as well which 
will be considered in Section 4 within a particular relevant example. 



2.1 The leptogenesis CP asymmetries 

We write the Lagrangian for the three left-handed lepton doublets Lq, the three RH 
charged leptons i^c, and the three heavy singlet Majorana neutrinos Ni as: 

- £ = \n,M,^N^ + L^Y^^NiH^ + L^Y^pEpHa (1) 

where H^^ are the SU{2) doublets Higgs fields, the matrix Mij is complex symmetric, 
and the Yukawa matrices Yai and Y^^ are arbitrary complex 3x3 matrices. 

The leptogenesis CP asymmetries for A^j-decays into a-leptons (a = e, /i, r) are defined 

as 

1 _ F. 



^In the cases when the flavor symmetry is broken close or below the leptogenesis scale, interesting 
scenarios that are quite difi^erent from the standard one, and that in general involve the dynamics of new 
scalar or fermion particles, might be realized [9]. 



where Fj — a and Fj — a denote respectively the decay rates into a-leptons and a-anti- 
leptons. We have defined Fj = Fj^Q, (with an analogous definition for Fj) so that the 
total A^j-decay widths are given by Fj + f j, and the total CP asymmetries ei = X^a^i^a 
are 

r. - r. 

(3) 



r, + F, 

Let us now introduce the Hermitian matrix 



(4) 



A one-loop perturbative calculation of the fiavor and of the total CP asymmetries e., 
and Ei gives [TO] : 

1 
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fij + Im 
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8ny 



(5) 
(6) 



where the hat denotes the basis in which the matrices M and of eq. (pO) are diagonal, 
with real and non-negative entries. The fij, Qi/s are known [10] functions of the heavy 
Majorana neutrino masses Mj. Indicating with Mi the lightest heavy neutrino mass 
{Ml = min{Mi, M2, M3}), defining Xi = Mf/Mf and introducing the function 
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2-x 



'1 + x) In 
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l + x 



X 



(7) 



such that ^(00) = —1, they are given by 



^ 3 i{xj/xi) 
2 ^Xj/xi 



and 



9ij 



In the limit of strongly hierarchical singlet neutrinos Xj ^ Xj, one has \gij\ <^ \fij\ ~ 
{Z/2)\J Xi/xj) and, barring strong phase cancelations in the Yukawa couplings, the second 
term in St^a can be neglected. In the quasi-degenerate limit Xj/xi ^ 1 we have \gij\ ~ 



\1 — Xj/ Xi\ ^ . It is apparent from eq. ([5]) that if 3^ is a diagonal matrix, all the 



fiavor dependent CP asymmetries vanish. As regards the total CP asymmetries in eq. ([6]) 
they vanish also in the case y has non-vanishing but real off-diagonal entries. 



2.2 Flavor symmetries 

We assume that the theory is invariant under transformations of a fiavor symmetry group 
G, spontaneously broken down to a subgroup if, through the vacuum expectation values 
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(VEVs) of a set of scalar fields (p. The ratios between {ip) and the ultraviolet cutoff A of 
the theory provide a set of small parameters rj ~ {(f) /A, in terms of which we can expand 
the quantities of interest. The symmetry group G can be discrete or continuous, global 
or local. The group G acts on the 3 left-handed electroweak doublets L, 3 RH neutrinos 
and 3 RH charged leptons E, as follows 

L' = nL{g)L , N' = nN{g)N , E' = nE{g)E , (9) 

where ^i{g) (/ = L, N, E) denote unitary representations of the group G for the generic 
group element g. After breaking of the group G, the Lagrangian eq. ([1]) is generated, and 
the three matrices Mij, Y^i and Y^^ will be in general functions of the symmetry breaking 
parameters 77. To study the implications of flavor symmetries for the leptogenesis GP 
asymmetries it is then convenient to expand these matrices in powers of t]: 

M = + 6M + ... , r = r° + + . . . , Y^ = Y^° + 5Y'' + ... (10) 

In these expansions M°, Y^ and Y^^ denote the terms surviving in the limit of exact 
symmetry 7] —>■ 0, while 6M, 6Y and 6Y'^ are all terms 0{ri). Flavor symmetries will imply 
different constraints on these matrices as well as on the matrix y of eq. (jlj), depending on 
their nature (non-Abelian or Abelian) and depending also on the specific representations 
to which the fields are assigned. We now analyze the different possibilities. 

2.2.1 Irreducible representations ^7^r 

We start by considering the case where f^Ar is a three-dimensional irreducible representa- 
tion of a non-Abelian group G. In this case the implications of the flavor symmetry for 
leptogenesis are simple and completely general: that is, in the limit of exact symmetry 
?7 — > all the leptogenesis GP asymmetries vanish. 

Let us consider the matrix 3^ of eq. (jl]) in the limit rj 0. Due to the invariance 
under G, = Y^^Y^ obeys the relation: 

nUg)y'^Nig) = y' (11) 

for any group element g. Therefore y^ commutes with ^lN{g) for any g. Then, by the 
first Shur's lemma, 3^*^ is a multiple of the identity matrix 3^" = \y\'^I, and Yq itself is 
proportional to a unitary matrix. This property is clearly basis-independent and holds in 
particular in the hatted basis, leading to the vanishing of all the GP asymmetries in the 
limit of exact symmetry: 

Ei^a = £i = for 7] = 0. (12) 
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Clearly, to accommodate the three A^'s in an irreducible representation, G should be non- 
Abelian, since Abelian groups have only one-dimensional irreducible representations. To 
study the dependence of the CP asymmetries on t] when the flavor symmetry is not exact, 
we expand the matrix y in the hatted basis obtaining: 

% = \y\%j + = y'^^Y + SY^Y'' + ■■■ (13) 

Inserting this in the eq. ([5]), and assuming for the time being that at the leading order fji 
and gji in eq. ([H]) are unrelated to r], we obtain: 



Im 



[Y:XjSy^j\ -0{v). (14) 



This is because while 6y is 0{ri), Y*^Yaj is not related to r] when no sum over a is taken. 
In contrast, by summing over a we obtain 

e.-Yl^m[6yf^ ^Oiv'). (15) 

Thus we can conclude that when the heavy Majorana neutrinos Ni belong to an irreducible 
representation of a non-Abelian flavor symmetry, the following two possibilities can be 
realized: 



1. If Mi < 10^2 dH [121 [131 [H [15] then the role of lepton flavor dynamics 
cannot be neglected in a description of leptogenesis from Ni-decays and the final 
asymmetry is a linear combination of the flavored CP asymmetries and therefore of 
first order in the symmetry breaking parameter t]. 

2. If Mi > 10^^ GeV (see footnote), then the one flavor approximation is accurate, the 
relevant quantities are the total CP asymmetries Ei and therefore the final asym- 
metry gets suppressed by one additional power of r] with respect to the previous 
case. 

In the cases when the non Abelian symmetry also implies degeneracies for the singlet 
neutrinos mass eigenvalues (|Mj — Mj\ -C Mj + M,), then we have gij ~ fij ~ 0{r]^^) and 
eqs. ([T^ and (ITBI) become Ei^a ~ C(?7°) and Ei ~ 0{ri). 

The previous analysis applies to several popular models that predict TBM, based for 
example on the discrete symmetry groups A4 pn[T8] and Zj x [19] (that were analyzed 
in [8] in the context of leptogenesis), T' [20], 5*4 [21], -D(4) [22], on the continuous group 
0(3) that occurs e.g. in models of minimal flavor violation [23], etc. 



^Notice that this condition becomes, more generaUy, Mi ^ 10^^ (l+tan^ f3) GeV in the supersymmetric 
case [TBI. 



2.2.2 Reducible representations Q 



N 



In the cases when the heavy Majorana neutrinos Ni transform according to a reducible 
representation Qn of the group G, it is still possible to derive some general conclusions. 
When at most one eigenvalue of M° vanishes, 3^" can be brought into diagonal form if 

yO^M^ - M^y° = 0. (16) 

Clearly, this implies that also in this case the leptogenesis CP asymmetries vanish in the 
limit of exact symmetry. Indeed, if 3^" is diagonal in the basis where M° is diagonal, 
then there exists a unitary matrix U that diagonalizes both matrices: M° = U*M^U'^ 
and = Uy^U'^ . Then, from [M°,3^''] = 0, the condition (fTGll immediately follows. 
Conversely, let us assume that eq. (fT6l) is satisfied. Since M° (like M) is a symmetric 
matrix, there exist (Takagi factorization) a unitary matrix U such that 

M^ = U*M^U^, (17) 

where M° is diagonal with real and non-negative eigenvalues. By acting with and 
U respectively on the left and right-hand sides of eq. ( fT6ll (and recalling the definition 
y^ = uy°W) we obtain 

(3>°)^M° - M°3>° (18) 
Since y^ = y*, eq. (ITS!) implies for the (ij) component of the matrix product above 

Myi) M - ^S] = 0' ^^(yi) M + = 0- (19) 

From eq. (fT9l) we can conclude the following: 

1. If the eigenvalues of M° are all different then 3^*^ is diagonal. For the leptogenesis 
CP asymmetries, the same conclusions stated at the end of the previous section 
hold. 

2. If two eigenvalues of M° are equal and nonvanishing, then 3^*^ can always be rotated 
into diagonal form. To prove this, let us assume e.g. that M^^ = Mj ^ 0. Then 
3^*^ is block diagonal with real entries in the corresponding 2x2 block 3^('jj)lj It 
is then possible to diagonalize y'^ij-^ by means of a 2 x 2 orthogonal transformation 

= -^2x2 while leaving M^^.^ = Mj /2x2 unaffected. 



■^Note that while this would imply that the total CP asymmetry Si vanishes, it would not necessarily 
imply the same for the flavor asymmetries Si^a, see eq. (1141) . 
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3. If all the three eigenvalues of M° are equal and non- vanishing, then, following the 
argument outlined above, can be always brought into diagonal form. 

Since in cases 2. and 3. the eigenvalues of M are characterized by degeneracies when the 
limit 1] is taken, the conclusion that the off diagonal terms of are at most of 0{ri) 
by itself is not sufficient to infer an estimate of the size of the CP asymmetries. As was 
already mentioned above, when Mj = Mj + 0{ri) we have that fij, Qij ~ 0{r]^^). Thus, 
only the total CP asymmetry Si gets parametrically suppressed by one power of t], while 
for the flavored CP asymmetries there is no parametric suppression. 

The last possibility, that is not contemplated in 1.-3., is that two eigenvalues of 
(or all the three) vanish, implying that two (three) light neutrinos are massless in the 
symmetric limit. In this case we have no information on the corresponding block y^^j-^ 
(on y^), and to reach some conclusion, we need to analyze the structure of the correction 
SM. For example, if M° itself vanishes, and 

y^^ 6M -6My^ = 0, (20) 

then again we can conclude that the off-diagonal terms 3^°^^- are at least of 0{ri). It is 
then clear that, in several relevant cases, in order to derive general conclusions about the 
size of the leptogenesis CP asymmetries the criterion eq. ( fT6l) , can still be applied to terms 
of higher order in t] in the expansions of M and 3^. 

Popular models for which the conclusions of this section can be relevant are based on 
the discrete symmetries S3 [2lj, -D4 [22], A* — ^ permutation symmetry. 

2.2.3 The Abelian case 

A particular case of a reducible representation Q^lg), is when it corresponds to the sum 
of three singlet representations: 



/ ioiig) \ 



(21) 



UJ2{9) 
V ujsig) J 

where [^^^(^f)^ = 1. If this occurs, in practice, the action of G on is Abelian. 

We can examine all possible cases. For simplicity, we assume that the transformation 
properties under G of the left-handed doublets are such that the maximum allowed 
number of entries in Y are different from zero. 

1. uji{g) = l (2 = 1,2,3). 

All the elements of M° are non-vanishing and are free-parameters of the model. 
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and so are the eigenvalues matrix M". In general, also the elements yi^j will be 
non-zero, since the matrix elements of M° and those of 3^ are completely unrelated. 
Concerning the CP asymmetries, the only conclusions that can can be inferred is 
that they will be generically non-vanishing. Models based on Abelian symmetries 
of the Froggatt-Nielsen type often realize this situation. 

2. uj2{g) = uJsig) = 1 and uji{g) =^ 1. 

is block diagonal with non- vanishing free parameters filling up M^^ , and con- 
dition eq. f|T6l) in general is not realized. If ui = —1, then M^^ ^ 0, otherwise 
vanishes. The invariance condition VL^y^^^t^ = y^ implies yf^ = 3^?3 = 0, 
and therefore all the CP asymmetries for A^^i vanish. The A^2,3 CP asymmetries are 
generically different from zero. 

3. ujsig) = 1 and uji{g), uj2{,g) ^ 1- 

M° is block diagonal with M^^ 7^ 0. If uoi = UO2 = —1 the M^^^ 2) block is non- 
vanishing with free parameters filling up all the entries. Clearly (modulo the 
replacement Ni N^) this case reduces to the previous one, and in particular 
Vt^j^y^VtM = y'^ implies y^i = y^2 — O vanishing of the N3 CP asymmetries. 

If ui, UJ2 7^ ±1 but co'iCo'2 = 1, then the 2) block is anti-diagonal. These two 
conditions also imply ujIuj2 7^ 1 and thus in the flavor basis ^ is a diagonal matrix. 
However, since condition (fT6!) in general is not realized, in the hatted basis y^ 7^ 
and the CP asymmetries for do not vanish. 

4. If uji 7^ ±1 (for i = 1,2, 3) and UiUj 7^ 1 for all i 7^ j, then = and the heavy 
neutrino masses arise at 0{6M) or at higher order. Unless uj*ujj = 1 for some i 7^ j, 
y is diagonal in the flavor basis, but condition ( JT8l) is generally not satisfied. Thus in 
the hatted basis 3^^^^- 7^ and all the CP asymmetries are generically non- vanishing. 

3 The flavor symmetry A4 

In recent years, the non-Abelian discrete group [17] stemmed out as one of the most 
interesting candidates for explaining the measured values of neutrino mass square differ- 
ences and mixing angles within a TBM [7] phenomenological pattern. In models based on 
A4 the heavy Majorana neutrinos are generally assigned to an irreducible representation 
of the group. Thus the main conclusion of Sect. 12.2.11 that is that Si, Si^a vanish in the 
limit of exact symmetry, applies [8] . For this reason A/^ represents a relevant case in which 
the constraints on the leptogenesis CP asymmetries implied by the flavor symmetry are 
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highly non-trivial. As we will discuss in this section, besides this general property, the 
specific realization of A4 on which we will focus, and that is based on the flavor sym- 
metry A4 X Z3 X U{1)fn, provides a highly constrained framework that results in some 
remarkably precise predictions for the neutrino sector. 



3.1 The discrete group A4 



We start by reviewing briefly the properties of the A4 group. A4 is the alternating group 
of order 4, e.g. the group of the even permutations of four objects, so it has 12 elements. 
Prom a geometrical point of view, it is the subgroup of the three-dimensional rotation 
group leaving invariant a regular tetrahedron. All the elements of the group can be 
expressed in terms of only two elements of the group itself (the generators) which we will 
call S and T. They obey the following rules: 



{STf 



(22) 



A4 has four inequivalent irreducible representations (irreps from now on) : three of dimen- 
sion 1 (1, V and 1") and one of dimension 3 (3). 
The form of the generators in the different irreps is given by: 



1 

1' 
1" 



3:^=3 



/-l 
2 

V 2 



S=l, 
S = l, 

2 2 \ 
-1 2 
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T 
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(23) 
(24) 
(25) 

(26) 



The product of two triplets decomposes as follows 

3 X 3 = 1 + 1' + 1" + 35 + 3a 
Explicitly, if we have two triplets a = (ai, 02, as) and h — (61, 62, ^3), the product reads 



{ah)k = ^ttiA^ijhj 



(27) 



^,3 



where k — 1, 1', 1", 35, 3^. The A^ matrices are given by: 





/ 1 









( ° 


1 


0^ 




/o 





1 ^ 




A^^ 












1 













1 





(28) 






1 











1^ 









0^ 





11 



1 
3 



/ 2 


^ 


a3s _ 1 


( ° 


-1 0] 




( ° 





-1 \ 





-1 


-1 





' 2 3 





2 





Vo -1 


) 




V 


2) 




v-1 





/ 



(29) 
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(30) 



where in eqs. f l29|) and fl30|) we have used the convention that the subscript labels the 
component of the triplet, e.g. 



(ab) 



3a 



]-{aAl^b, aAl'^b, 0^3^ 6) = ■^(a2&3 - 0362, ^1^2 - a2&i, cts^i - ^ih 



(31) 



The group has two obvious subgroups: Gs — Z2, the reflection subgroup generated by 
S, and Gt — Z^, which is generated by T. It is immediate to see that the VEVs 

(¥.t)cx (1,0,0) 
{<fs)^ (1,1,1) 

break A4 respectively to Gt and Gs- As we will see in the next section, both triplets v^t 
and (ps are needed in order to get the TBM matrix for the neutrinos. 



3.2 The A4X Z3 X U (1) FN model 

We recall here the main features of a well known model [6] for TBM that is based on 
the symmetry group x Z3 x ?7(l)_FAf B In order to write the most general Lagrangian 
invariant under the symmetry group, we first have to assign the fields to specific A4 irreps. 
We assume that the model is supersymmetric, this turns out to be useful to ensure the 
correct vacuum alignment of the scalar fields breaking A4, that is a crucial feature to 
obtain TBM. We divide the field content of the model in three different sectors: 

1. matter fields sector: the lepton doublets La transform as a triplet 3, while the RH 
charged leptons Ei = e'^, E2 = fi'^ and = t'^ transform respectively as 1, 1" 
and 1'. In order to discuss leptogenesis, we need to introduce also the three heavy 
Majorana SU{2) singlet neutrinos Ni, which are assigned to another triplet 3. 

^Several models based on the A4 symmetry have been constructed. Recently the baryon asymmetry 
was computed in one of these variants, where the RH neutrino masses are strongly degenerate |25) . 
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2. symmetry breaking sector: the Higgs fields Hu and Hd (with opposite hypercharge) 
transform in the 1 irrep of A4. We introduce four fiavons for the spontaneous 
breaking of A4 (pT — "i, — "i-, i — '^ and ^ ~ 1. 

3. driving sector: we introduce three fields, v?o 5 V^o ' ^O; that allow to build a non trivial 
superpotential that, once minimized, will ensure the correct vacuum alignment. 

On top of we also impose a discrete symmetry, that is introduced in order to guar- 
antee that, at the leading order, the charged leptons and the neutrino sector selectively 
couple to two different sets of flavons: ipx for charged leptons and v^s,^,^ for neutrinos. 
This ensures that the vacuum alignment produces the tribimaximal mixing matrix for 
the neutrinos [6]. The symmetry group also includes a U{1)fn Froggatt-Nielsen factor, 
that acts only on the RH charged leptons. It is spontaneously broken by the VEV of 
a flavon field 9 carrying a negative unit of the Froggatt-Nielsen charge. By choosing 
U{l)p]^ charges (2,1,0) for (e'^, /i"^, r^), the mass hierarchy in the charged lepton sector 
can be reproduced. Besides the fiavor symmetry, a continuous U{1)r i?-symmetry, con- 
taining the usual R-parity of the SUSY models is also needed. Matter, symmetry-breaking 
and driving fields have respectively i?-charge 1, and 2. The fields assignments to the 
/I4 X Z3 X U{1)fn irreps are summarized in Tabled! 





L 










Hu,d 


iprp 
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T 








3 


1 


1" 


V 


3 


1 


3 


3 


1 


1 


1 


3 


3 


1 


^3 








a;2 




1 


1 








1 


1 






U{1)fn 





2 


1 























-1 











U{1)r 


1 


1 


1 


1 


1 




















2 


2 


2 



Table 1: The fields of the A4 x Z3 x U{1)fn naodel and their representations. 



The most general superpotential compatible with the representation assignment of 
Table [T] is given by 

W = Wl + Wd 

wl is the leptonic part of the superpotential, 

WL = ye[^L)^e'Hd + y^.[^L)^^i'^Hd + yr[^L)^T-Hd 

+ y{LN)^Hu+(xa^ + £a(){NN)^ + x,{^sNN)^ + ... (32) 

where the dots stand for higher dimensions operators, suppressed by additional powers of 
the cutoff A. We note that it is always possible to redefine the fields in a way that ?/e, y^, 
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Ur and y are real numbers, while Xq, Xa and Xh are in general complex parameters, is 
the "driving" part of the superpotential 

+ 91 i^o'PsVs) 1 + 92^ iVoVs) 1 + ^3-^0 {^sVs)i 

+ 94^0^'' +95^0^^ +96^0^'' + - 

where dots denote higher dimensional contributions. For a detailed discussion about the 
minimization of the superpotential we refer to [6]; here we only mention that, given the 
previous expression for wd, the VEVs of the fiavon fields at the leading order are 



(0 



vt{1,0,0) 

^^5(1, 1,1) 





(33) 



where vt, vs can be expressed in term of m, g and Qi {i = 1 ... 6) while u remains 
undetermined. We assume that all VEVs in the fiavon sector are of the same order: 
Vt ^ Vs ^ u ^ V. Their value in units of the cutoff scale A, V/A, corresponds to the 
parameter t] of the previous section §. 

From eq. ( l32l) we see that already in the phase when the flavor symmetry is unbroken, 
the Yukawa matrix for the singlet neutrinos is non-vanishing: 



/ 1 \ 
1 

Vo 1 o; 



(34) 



In contrast, as long a.s A4 x Z3 x U{1) fn is, unbroken, the Yukawa matrix for the charged 
leptons and the Majorana neutrinos mass matrix both vanish. In terms of the expansions 
eq. (|TOl) we thus have = M° = 0. After the breaking of the flavor symmetry we get 



f Ve 



V 



\ 









Vt ) 



M = 6M 



I a + 26/3 

-6/3 
V -V3 



-6/3 
26/3 
^-6/3 



Vt 

X' 



-6/3 \ 
^-6/3 
26/3 



(35) 



(36) 



^The precise definition of 77 in this model will be given in section 4. 



14 



where a = 2xaU and b = 2xbVs are complex numbers. The complex symmetric matrix M 
is diagonalized by the transformation 



M = UMU^ 

where M is a diagonal with real and positive entries given by 

M = diag(Mi,M2,M3) = diag(|a + 6|, \a\, \-a + b\). 
The unitary matrix U can be written as 



where Utb is the TBM matrix 



TB 



TB'^PH, 



V5 7! ° 

1 1 1 

1 1 1 

\ Ve Vs V2 J 



(37) 



(38) 



(39) 



(40) 



and UpH = diag(e*"i/^, e*"2/2^ e*a3/2^ jg ^ matrix of phases, with ai = — arg(a + b), 
a2 = — arg(a) and as = — arg(— a + b). 

After electroweak symmetry breaking {Hu,d) = Vu,d the mass matrix for the light 
neutrinos is obtained from the seesaw formula: 

= -vlY'^M-^Y'''' = -vly^M-K (41) 

Since from eq. (|371) U*M^^U'^ = diag(Mf^, Mg"^, M^'^, ), we obtain for the light neutrinos 
mass eigenvalues 

mi = ^, (42) 

where we defined 



yp = y sin (3 



tan/5 



Vd 



(43) 



and V = \/v\^v\ ~ 174 GeV. Notice that these mass relations correspond to a particular 
model of sequential dominance [26]. The neutrino mixing matrix is given by 

U, = U'^ = UtbUpu ■ (44) 
Thus, at first order in the flavor symmetry breaking, the neutrino mixing matrix is TBM. 
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3.3 The neutrino mass spectrum 



According to eq. (H2ll . in the present approximation the hght neutrino mass spectrum is 
directly related to the heavy neutrino masses. These, in turn, can be expressed in terms 
of just three independent parameters (cf. eq. fl38l) ) that for example can be chosen to be 
\a\ = M2 = y'^v'^/m2, \z\ and ip. The latter two are defined according to 

- = |2|e'^. (45) 
a 

Experimentally, only two observables related to the spectrum have been measured. For 
normal (inverted) hierarchy they are [27j 

ml -ml ^ AmL = (7.67l°f ) X 10-^ eV^ , 
\ml-ml{ml)\ = Am^^^ = (2.46 (2.45) ± 0.15) x 10^^ eV^ , (46) 

and thus the neutrino mass spectrum is not fully determined. However, within the present 
model the value of a third parameter is bounded | cosy^l < 1. We now show that this 
condition constrains significatively the absolute neutrino mass scale. A relation between 
the phase ip and the neutrino masses is easily derived. We write 

w»=^ = 1(4-4). 



'1 / ml m^ 

2 \m\ 777,3 

and taking the ratio of the previous two equations we obtain 





to to 


_ ml 






2 

mi 


ml 






/ml 
\ml 


ml) 


- 1 



'"'1 "^■^ /^^\ 

COS(y9= ^ . (49) 



This result holds both for normal ordering (NO) and for inverted ordering (10). By 
expressing the heavier masses in eq. fH^ in terms of the lightest one 777^ = 7771 (7773) for 
NO (10), and of 777soi = \/ ^ral^i ^ 0.0088 eV, m^^m = \/ ^ml^^ ~ 0.050 eV and solving 
for the condition | [cosy?] (777;) | < 1, we obtain, in agreement with [28l [29], the following 
limits for the lightest neutrino mass: 

0.0044 eV< mi < 0.0060 eV (NO) (50) 
m-i > 0.017 eV (10). (51) 
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Nomal hierarchy Inverted hierarchy 




0.0044 0.0048 0.0052 0.0056 0.006 0.0167 0.02 0.025 0.03 0.035 



Figure 1: cosip versus the lightest neutrino mass for the NO (left) and 10 (right) cases. For 
the NO plot the most conservative limits are obtained for Am^^„ = (2.46—0.15) x 10~^ eV^ 
together with Am^^^ = (7.67 — 0.21) x 10~^ eV^ (left red curve), and with Am^^; = 
(7.67 + 0.22) X 10-^ eV^ (right blue curve). For the 10 plot the most conservative limit 
is obtained for Am^^^ = (2.45 - 0.15) x 10"=^ eV^ and Am^^^ = (7.67 - 0.21) x 10"^ eVl 



The function [cos !f]{mi) is plotted in Figure [T] for the NO (left) and 10 (right). By ex- 
panding [cos (p]{mi) in powers of r = mgoi/matm, we can also derive approximate analytical 
expressions for the limits on mi. For NO we obtain: 

-^i^l-^r + ... j < < (^1 + + ... 1 , (52) 

and for the 10: 



where the dots represent higher order terms in the expansion in powers of r. 

For NO we have both a lower and an upper bound on mi, that select a rather small 
range for the possible rrii values of width ~ ^/3r and centered around nisoi/VS. Thus, 
the neutrino mass spectrum is essentially determined. From Fig. [1] we can also see that 
the phase ip always remains quite close to zero (cosy^ < 7r/15), and that both the upper 
and lower bounds are saturated for ip = 0. In the 10 case, we only get a lower bound on 
7713, that is saturated for ip = ±7r. The neutrino mass remains unbounded from above, 
and the phase if is allowed to vary between 7r/2 and vr or between — vr and —tt/2, that is 
in the ranges where cos ip is negative. 

The results presented so far are of course approximate since the model gets corrections 
when higher dimensional operators are included in the Lagrangian. The inclusion of higher 
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Figure 2: Plots of the heavy (upper panels) and light (lower panels) neutrino masses 
divided by function of the lightest left-handed neutrino mass, according to eq. ( H2l) . 

In the left panels (NO case) the two vertical green lines show the allowed mi mass range 
as given in eq. (150|) . In the right panels (10 case) the vertical green line shows the lower 
limit on ms eq. (I5T1) . The dashed area is the region excluded at 95% CL by the upper 
bound (cf. Eq. (IFH) ) obtained by combining the cosmological upper bound on the sum of 
neutrino masses with neutrino oscillations data. 



dimensional operators has also the effect of shifting the VEVs of the flavon fields from their 
leading order values, eq. (133|) . By assuming that the VEVs of the flavon fields have similar 
values vt ^ Vs ^ u ^ V , these corrections modify the leading order approximation by 
terms of relative order V/ K. The allowed range of V/ K is determined by the requirement 
that sub-leading corrections which perturb the leading order result are not too large and 
by the requirement that the r Yukawa coupling Ut does not become too large. The first 
requirement results in an upper bound on V/ K of about 0.05, which mainly comes from 
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the fact that the solar mixing angle should remain in its la range. The second one gives 
a lower bound which we estimate as 

V ta.n 3 rrir tan/? , 
-—^0.01 . 54 

A yr V yr 

By asking \yr\ < 3 we find a lower limit on V/A close to the upper bound 0.05 for 
tan/3 = 15, whereas tan/3 = 2 gives as lower limit V/A> 0.007. We choose as maximal 
range: 

0.007 < ^ < 0.05 , (55) 

which shrinks when tan (3 is increased from 2 to 15. 

In particular the leading order expressions of neutrino masses are modified by terms 
of relative order V/A. However, close to cosip = ±1, where the bounds are saturated, 
the corrections to both the numerator and the denominator of eq. ( H9i) remain of relative 
order V/A, and thus the bounds in eq. fl5U]) and eq. S5T\i are not significantly affected. 
For normal hierarchy, that requires cos (f very close to one in the full allowed mass range 
of eq. (150|) . the leading expression for cosy? given in eq. (l49l) always remains a good 
approximation. In the case of inverted hierarchy, close to cos if = —1, when is near 
its lower bound, the corrections are also negligible. Deviations from eq. (H9l) can become 
significant when cosip approaches —0.2. This happens for ~ 0.09 eV and mi^2 ~ 0.10 
eV, that is when the spectrum becomes nearly degenerate. 

It is interesting to estimate the order of magnitude for the RH neutrinos masses. Since 
no suppression is expected for parameters that are unrelated to the breaking of the fiavor 
symmetry, we take = 0{1). With this, and using as the light mass scale in the seesaw 
equation (H2l) matm — 0.05eV, we obtain 

Mi ~ lO^^^^^GeV . (56) 

A detailed summary of the neutrino mass relations and bounds is given in figure [21 where 
we have plotted the three light and the three heavy neutrino masses, taking into account 
both the information from neutrino oscillations data (cf. eq. (H6ll ) and the seesaw relations 
eq. In this way, at leading order, all six masses can be expressed as a function of just 
one independent parameter, that can be conveniently chosen to be the lightest left-handed 
neutrino mass m^. Of course there is also a dependence on the neutrino mass ordering, 
namely if m; = mi or rrii = m^. In the figures we also display the cosmological upper 
bound on mi 

mi < 0.2 eV (95% CL) , (57) 
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that follows from the WMAP5 |30] upper bound on the sum of the three neutrino masses 
combined with the constraints from mass squared differences from oscillations data (cf. 
eq. 



3.4 Neutrinoless double-/? decay 

In the approximation of neglecting terms of higher order in the symmetry breaking, as 
well as RGE running effects from the high scale to the eV scale, we can straightforwardly 
obtain predictions for the Oz/2/5 decay parameter \mee\ for both the NO and 10 cases. The 
0u2P decay parameter is defined as 



me 



(5i 



and corresponds to the (11) entry in the neutrino mass matrix m,^ in eq. ( HTl) : 



\mee\ = = vly"^ 



+ 



1 



(59) 



3(a + h) 3a 

In terms of physical neutrino masses, and to the order we are working here {{Uu)id, = 0) 
eq. (J58l) assumes the particularly simple form: 



irrip 



3 3 



(60) 



and thus depends only the phase difference 0:2— ai, where a, are the phases of the diagonal 
matrix Uph defined below eq. fj^Ol) . 

In order to compute the allowed range of \mee\ when the lightest neutrino mass is 
allowed to vary in the allowed region eq. (150|) (for NO) or eq. ( l5Ti) (for 10), it is more 
convenient to express |mee| directly in terms of the phase cos0 in eq. (149|) . This can be 
done more easily by using eq. (159|1 . and yields: 



rrir. 




1 + 4 



2 + 


z 


COS0 


1 + 


z 


2 + 2|z 


COS (f) 



(61) 



that holds for both NO and 10. 

Following [21] we parameterize the forecast sensitivity of future 0i/2/3 experiments as 



10 h meV, 



(62) 



where h = 0.6 -r- 2.8 parameterizes the theoretical uncertainty related to different nuclear 
matrix elements calculations. We then confront the prediction of the A4 x x U{1)fn 
model with eq. (1^ . 
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Figure 3: \mee\ versus the lightest neutrino mass. Left paneh the NO case. The thick 
blue line depict the value of |mee| as a function of mi. The thin red line illustrates the size 
of the corrections that could arise for sin^ Oi^ ^ 0.026. Right panel: |mee| as a function of 
772.3 in the 10 case. In both cases the dark shaded horizontal region depicts the projected 
sensitivity of future Oz/2/3 experiments. 



For the NO case we have m; = mi and, according to eq. (I3U]) . m; can vary between 
0.0044 eV and 0.0060 eV. The results for |mee| in this case are depicted with the thick 
blue line in the left panel of figure [31 The dark horizontal region corresponds to the 
expected sensitivity of future 0^2^ experiments according to eq. fl62l) . In the leading 
order approximation only mi and m2 contribute to |mee| (see eq. ( 160|) ). In the NO case 
these are the two smaller masses, and it is then reasonable to ask to what extent effects 
related to the largest mass ms, that also contributes to |mee| when sin^ ^13 7^ 0, can 
affect these results. By inserting in Uu the maximum value sin^ ^13 ^ 0.026 allowed at 
la by present data [27], with the value of m3 fixed in terms of mi and the mass squared 
differences, we obtain the values of |mee| depicted in figure [3] (left panel) with the thin 
red hne. Of course, the procedure of setting {Uu)i3 7^ by hand, and not as a result 
of the inclusion of higher order terms, can just give a feeling of the possible effects of a 
non- vanishing sin^^^i3, but does not correspond in any way to an improved prediction, 
that would require a consistent treatment of all higher order effects. 

From the left panel in figure [3] we can then conclude that for all the allowed values 
of mi, even in the most optimistic situation, |mee| remains at best marginally in the 
reach of the future experiments sensitivity. We also see that higher order contributions 
oc m3 ■ sin^ ^13 are not expected to change this conclusion. 

For the 10 we have mi = m^. A plot of \mee\ as a function of m3 > 0.017eV 
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(see eq. (15T!) ) is depicted with a thick blue hne in the right panel of figure [31 We see 
that in the 10 case, |mee| remains well above the sensitivity of future experiments, ex- 
cept in the most pessimistic situation {h = 2.8) and when the lower bound is saturated 
(ms ~ 0.017 eV). In this situation the interference of the two contributions in eq. (l60ll 
is maximally destructive, yielding \mee\ — The results plotted in figure [3] (right 

panel) will not be significatively affected by contributions oc rris ■ sin^ ^13, that for 10 are 
surely negligible, or by other higher order effects. Thus, the A/j^ x x U{1)fn model 
predicts that if the neutrino ordering is inverted, a Oz/2/5 decay signal quite likely will be 
observed in the next future. Detection of Oz/2/3 decay is guaranteed if new theoretical 
computations will establish that h < 1.7. 

4 Leptogenesis in the A4 x x U{1)fn model 

In section [2.2. II we have shown that an important consequence of non-Abelian flavor sym- 
metries with the Ni assigned to an irreducible representation of the symmetry group, is 
that the leptogenesis CP asymmetries vanish at leading order. Their size is thus deter- 
mined by the size of the flavor symmetry breaking parameters rj. In the A^x Z^xU {1)fn 
the heavy Majorana masses are bounded from below by Mi/y'^ > 1.5 x lO^'' GeV (see 
upper right panel in figure El), and thus are generally quite large. Therefore, the condition 
Mi > 10^^ (1+tan^ P) GeV, that ensures that leptogenesis occurs in the unfiavored regime, 
is generally matched by all the for natural values of the parameters (e.g. tan/5 ^ 10 
and y > 0.5 (NO) or y > 0.8 (10)). In this regime, the relevant quantities that have to 
be considered are the total CP asymmetries e,, and since these asymmetries are 0{7f) 
(see eq. (1151) ). the requirement of successful leptogenesis can provide hints on the mini- 
mum size of the typical symmetry breaking effects. The important point here is that in 
general the symmetry breaking parameters are also related to some observables that are 
measurable at low energy. For example, in models that predict neutrinos tribimaximal 
mixing (TBM) (as is the case of our A^x Z^x U{1)fn model) both the value of ^13 and 
the deviation from maximal angle of ^23 are related to r] and, by studying leptogenesis, 
one can infer preferred ranges of values for these parameters. This is a new type of low- 
energy/high-energy connection (see however [8l[32]) and constitutes the main motivation 
of the following analysis. 
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4.1 The leptogenesis CP asymmetries at subleading order 

Since in our A4 model the heavy Majorana neutrinos belong to an irreducible representa- 
tion of the flavor group, independently of the particular basis we have that = y°Tyo _ 
y"^ I. Thus, at leading order, all the CP asymmetries vanish. When the leading symmetry 
breaking terms are introduced in the Lagrangian, new complex parameters are generated 
that give rise to non- vanishing CP asymmetries (see eqs. f[T^ and f[T^ ). 

All the possible leading order corrections to the ^4 x Z3 x U{1)fn model have been 
listed in [6]. However, as long as the CP asymmetry are concerned, the relevant terms are 
only two [8], and both represent corrections to the Majorana neutrinos Yukawa matrix: 



s) Hu-, 



a) Hu 



(63) 



Using the product rules of eqs. fl28H30l) . the Yukawa matrix now reads: 



Y = Y° + 6Y 



\ 






d 



y 





d 




(64) 



( y + 2d 


y 

where we have introduced the two complex numbers c, d (|c|, \d\ <^ y) defined as 

y2VT , yiVT 

c = — — , d = — — . 

2A ' 3A 

Assuming, as is reasonable to do, that yi and 1/2 are numbers of the same order, we can 
parameterize the size of the symmetry breaking effects with a single hierarchical parameter 
7], defining for example 

Re(c) 



(65) 



1] = Re((i) and 



P 



Re{d) 



(66) 



where p = 0{1). Then in the hatted basis the matrix y = Y'^Y becomes 

y = UpijU^syuTBUpH 



y + 2ri 



2^2 r^e*"^! -pr]e'"''' \ 



y-2ri 



1032 



(67) 



Since 3^ is a physical observable (it can be measured in principle by measuring the CP 
asymmetries) it depends only on phase differences a^j = (a^ — aj)/2. 
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Recalling now the expressions of the CP asymmetries given in eq. ([6]) in terms of 3^ 
and of the functions fij (see eq. ([8])) we can write: 

£i = — 8 sin(2a2i)/i2 + ttP^ sin(2a3i)/i3 

„2 / g 

£2 = ^ I -8 sin (2a2i) /21 + 77 sin (2a32) /23 

4 „2 ... N . 8 .2 



^3 = ^ --p'sin(2a3i)/3i-7rp'sin(2a32)/32 . (68) 
87r \ 3 3 / 

Starting from 2 021 = — ai'g(fl) + arg(a + b) and 2a32 = arg(a) — arg(— a + b), one finds 
the following relations 

\z\ simp 



+ 2 \z\ cos(p + I^P 



sin(2a;2i) 
sin(2 032) 

A/i - z 

, , 2 1^;! sin 09 , , 

sin(2«3i) = - „ ' ,1 ,2 . 2 • (69) 
\z\^ -1)2 + 4 \z\^ sm^v9 







sin 


VI -2 




cos + 


\z\ 


2 



We recall that, after expressing the heavier neutrino masses in terms of mgoi, 'n^atm and 
mi, one has that sin(2a2i), sin(2a32) and sin(2a3i) in eq. (169|) are a function of mi only. 
Moreover, given that the expressions for the depend, through the functions fij, only on 
ratios of the heavy Majorana masses Mj, and that through the seesaw formula eq. fHTl) 
these ratios are directly related to the ratios of light neutrino masses mi, we can con- 
clude that the CP-asymmetries in eq. (I68l) depend only on mi, on the non-hierarchical 
parameter p, and on the parameter 77 that quantifies the flavor symmetry breaking effects. 



4.2 Analytic approximations 

In this Section we present simple analytical formulae to estimate the matter-antimatter 
asymmetry produced via leptogenesis within the x Z3 x f/(l)p'7v model. The presence of 
the heaviest RH neutrino makes possible for the CP asymmetry £2 of the next-to-lightest 
RH neutrino to be unsuppressed (cf. (168|) ) compared to the CP asymmetry of the lightest 
RH neutrino, even when the RH neutrino spectrum is strongly hierarchical [33]. Further- 
more, since in our case the heavy neutrino spectrum is only mildly hierarchical, even the 
CP asymmetry of the heaviest RH neutrino is not particularly suppressed. Therefore, all 
the three contributions to the final asymmetry have to be taken into account since in the 
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presence of compensating effects from reduced wash-outs the B — L asymmetry generated 
in the decays of the heaviest RH neutrinos could become comparable to the one generated 
in the decays of the two lighter RH neutrinos. This situation is indeed realized in the NO 
case. 

For our estimates, we adopt the following simplifications: 

• For the reasons explained in the beginning of the Section, we assume that leptoge- 
nesis occurs in the unflavored regime. 

• The Z^^U{1)fn model must be supersymmetric, since the vacuum alignment 
conditions eq. ([3SD are fulfilled within a supersymmetric framework [6J . Our results 
are instead obtained neglecting all supersymmetric partners effects. This underes- 
timates the resulting B — L asymmetry by a factor ^ (see ref. [3l] and [15] Sec. 
10). 

• The value of the final B — L asymmetry that we estimate is obtained by summing 
up the asymmetries generated in the decays of the three heavy neutrinos [HI [37] 
but neglecting the wash-out of the asymmetry due to the inverse processes of the 
lighter RH neutrinos [2SIES]- This can be done because, neglecting 0(77^) terms, 
the leptons produced in the decays of the three RH neutrinos are orthogonal to each 
other: 

mL,)\' = = 6,, + 0{v') . (70) 
I JjjJ^mI 

Doing this we are slightly overestimating the final asymmetry since we neglect the 
fact that the three lepton states = (yu)'^^'^ X^a^ail-^a) produced in de- 
cays are not exactly orthogonal one to the other, and thus part of the asymmetry 
produced by the heavier states gets washed out by the lighter RH neutrinos interac- 
tions. These effects could be taken into account in a straightforward way following 
the procedure explained in in [33]. However, neglecting the wash-out from lighter 
RH neutrinos is certainly consistent with the order of our approximation since, e.g. 
relative 0{ri) corrections to the CP asymmetries, that could produce even larger 
effects, are also neglected. With this approximation we also do not have to worry 
about complications coming from an overlap between decays and inverse decays that 
occur when the RH neutrino mass spectrum is not strongly hierarchical, and the 
leptons produced in the decays of the three RH neutrinos are not orthogonal to each 
other [38] . 
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• We neglect subleading leptogenesis effects like AL = 1 scatterings [39l HQl HI HI] and 
CP violation in scatterings [l2], thermal corrections [4j, spectator processes [i3lHi] . 
departure from kinetic equilibrium [15]. In the strong wash-out regime that is the 
relevant one for our model, these effects give corrections at most at the level of 
~ 50%. 

With these approximations, the B ~ L asymmetry can be estimated by solving the 
following three independent pairs of Boltzmann equations 



^ = -A (iV^. - AT^^) , (^ = 1,2,3) (71) 



dN 



(0 

— r- r> fAT.. Ar(*) 



dzi 



A [Nn, - K^) - Nli^ [Wi{zi) + ATy,(z,)] , (72) 



where Zi = Mi/T. We indicated with Nx any particle number or asymmetry X calcu- 
lated in a portion of co-moving volume containing one heavy neutrino in ultra-relativistic 
thermal equilibrium, so that N'^,{T ^ Mj) = 1. With this convention the predicted 
baryon-to-photon ratio rjB is related to the final value of the B — L asymmetry Nb-l by 
the relation ^ 

VB = «sph%^ ^ 0.96 X 10-' Nb-l , (73) 

7 

where N^'^ ~ 37, and a^ph = 28/79. The decay factors are given by 

D,^^ = K,zA-\ . (74) 



Hzi \% 

Moreover, indicating with g^^ = qsm = 106.75 the total number of degrees of freedom and 
with Mpi = 1.22 X 10^^ GeV the Planck mass, the expansion rate can be expressed as 



mz-)- J^^--166 /V^- (75) 

The total decay rates, j = (Fj + Fj) (l/7i), are the product of the decay widths times 
the thermally averaged dilation factors that can be expressed in terms of the ratio of the 
modified Bessel functions, such that (l/7t) = /Ci(-2j)//C2(-2j) • The equilibrium abundance 
and its rate can be also expressed in terms of the modified Bessel functions: 



K^M) = \ ^\ ^2{z,) , ^ = ^^(^^) • (76) 
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Introducing the effective washout parameters 



mi = V — — 77 



and the equihbrium neutrino mass [l6l [3] 

^Ts Mpi 
the decay parameters can be expressed as 



= '"^ ~ 1.08 X 10"^ eV , (78) 



= 1) m* 

In our case, from the Eq. ( !67l) one can verify that mj ~ mj and therefore there is a very 
simple relation between neutrino masses and decay parameters. 

After proper subtraction of the resonant contribution from AL = 2 processes [17], the 
inverse decay washout terms are given by 

W,{z,) = ^K,IC^{zi)zf. (80) 

The wash-out term AWj(zj) is the non- resonant contribution to the wash-out coming from 
AL = 2 processes and can be written as 



AW,{z,)--M,ml, (81) 



where 

_ 3v/5Mpi 
"-4C(3)vr9/2t;4y^- ^''^ 

The B — L asymmetry produced from A^j-decays can then be estimated as [5] 

Nf_^= e,K{K{)e-^^''''^l (83) 

where K{Ki) accounts for the wash-out from inverse processes and is approximately given 
by 

2 / KiZg{KA\ 



The quantity 



^B(ir,) ~2 + 4iri°-i=^e-^ (85) 



gives the approximate the value of Zi around which the final asymmetry from Aj-decays 
is dominant ly produced. 
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The exponential factor in eq. (183|) accounts for the wash-out from AL = 2 processes. 
One can notice that the two wash-out contributions factorize. Notice also that the AL = 2 
processes suppression is relevant only for Mi > 10^'^GeV (0.05 eV/mj)^. 

Finally, using eq. ( |83ll and eq. (1731) . the baryon to photon number ratio can be written 

as 

= ^riic^ 0.96 x 10"^ ^ eiK{Ki) e~^^^'"'\ (86) 

i i 

This theoretical prediction has to be compared with the observed value from WMAP data 

m 

r/^^^ = (6.2 ±0.15) X 10-1°. (87) 

4.3 Results 

We now describe the results for 77^, separating the discussion for the NO and the 10 case. 
4.3.1 Normal Ordering 

In the upper panels of Fig. H] we show the dependence on mi = mi of the three CP 
asymmetries Ei divided by the square of symmetry breaking parameter rj'^ for positive 
values of sin (p. For negative values they are simply all opposite. Therefore, by switching 
the sign of sinyj the sign of the final baryon asymmetry changes as well (i.e., the model 
does not predict the sign of the baryon asymmetry). It turns out that the correct sign of 
the baryon asymmetry is obtained for sin (p < 0. In the lower panels of Fig. H] we show 
the absolute values of rji/r]'^ {i = 1,2,3), that is the relevant quantity, together with the 
total value |?7b| = | Y^^fjil- We show the plots for three different values of p = 0.5,1,2 
from left to right and for three different values of = 1, 2, 3 from above to below. One 
can see how for increasing p typically the finally asymmetry increases. 

On the other hand, increasing yp, the three RH neutrino masses increase up to a 
critical value above which there is an exponential suppression from AL = 2 processes. 
This critical value is about ~ 2 and for = 3 the value of rjB/rj^ reproduces the 
observed asymmetry for the highest possible values of the //-range 5 x 10^^ -^ 5 x 10~^. 

It is also interesting to notice that the contribution from the lightest RH neutrinos 
(depicted in Fig. H] with the green dotted line) is sub-dominant, despite the fact that the 
lightest RH neutrino CP asymmetry l^al is typically the highest or anyway comparable 
to \e2\ (depending on the value of mi). The reason is that ~ matm/"^* ^ -^1,2 and 
therefore the wash-out is much stronger compared to the two heavier RH neutrinos. 
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Figure 4: NO. In the 3 upper panels we plot the CP asymmetries Ei divided by r/^ as 
a function of mi for simp > (note the different vertical scale in the last panel). The 
dashed line corresponds to i = 1, the short-dashed line to z = 2 and the dotted line to 
i = 3. The solid line is 1 — cos(f (cf. Eq. ( l49l) ). In the 9 lower panels we have plotted the 
three \'r]i\/f]'^ (cf. Eq. ( l86l) ). with the same line convention as for the CP asymmetries for 
i = 1,2,3 respectively, and \rjB\/if' (solid line). All the left panels correspond to p = 0.5, 
all the central panels to p = 1 and all the right panels to p = 2. The three rows of nine 
lower panels, showing the asymmetries, correspond to yp = 1, = 2 and yp = 3 starting 
from above. 29 
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Figure 5: NO. Plots of the values of i] that reproduce the observed r^^, as a function of mi 
for p = 0.5 (left panel), p = 1 (central panel), p = 2 (right panel). The different curves 
correspond to = 0.5 (dotted), ?//3 = 1 (short-dashed), yp = 2 (long-dashed), yp = ?> 
(dot-dashed) and yp = A (solid). The gray band is the indicative optimal range of values 
of = 5 X 10-3^5 X 10-2. 



In Fig. Owe have plotted, as a function of mi and for different choices of p and yp, the 
value of 7] such that rjB = t]'^^^^. One can see that this value always falls in the optimal 
range r] = 5 x 10^^ ^ 5 x 10^^ (the grey band). It seems therefore that, despite the many 
constraints on the model parameters and in particular the fact that there is only one 
independent complex phase and that the single dimensional parameter mi is practically 
fixed, the model reproduces the observed baryon asymmetry for natural values of the 
parameters in quite a satisfactory way. 

4.3.2 Inverted Ordering 

In the upper panels of Fig. [6] we show the dependence on = mi of the three Ei/r]'^ 
for positive values of sin if. Again, by switching the sign of sin if, one can always change 
the final sign of the final asymmetry and again the correct (positive) sign is obtained 
for sin If < 0. We show again three different values of p = 0.5,1,2 from left to right. 
However, one can see that now there is only a tiny dependence on p. The reason can be 
easily understood inspecting the three Eq.'s (168|) . The strong degeneracy between the two 
lower RH neutrino masses. Mi and M2, implies |^(a;2/a;i)|, |^(a;i/a;2)| ~ 10^ and therefore 
1/12 1; 1/21 1 ^ 1- In this way the term depending on p gives a negligible contribution. 

We again show examples for three different values of y^j. This time we choose y^ = 
0.5, 1, 1.7 from above to below. One can see how for = 0.5 one has rjB/rf ~ 10~^-^10-'^ 
and therefore too small values of < 10"'^ are required to explain the observed asymmetry. 
On the other hand for = 1.7 one has rjB /iff ~ 10~^ -r- 10~^ and this time too large 
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Figure 6: 10. In the 3 upper panels we plot the three Si/rf as a function of mi and for 
sin 99 > 0. The dashed line corresponds to z = 1, the short-dashed line to z = 2 and the 
dotted line to z = 3. The solid line is 1 + cosip (cf. Eq. fH9l) ). In the 9 lower panels we 
have plotted the three \'rii\/rf (cf. Eq. ( 1861) ) with the same line correspondence as for the 
CP asymmetries and \riB\/ff (solid line). All the left panels correspond to p = 0.5, all the 
central panels to p = 1 and all the right panels to p = 2. The three rows of nine lower 
panels showing the asymmetries correspond to yp = 0.5, yp = 1 and yp = 1.7 starting 
from above. 
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Figure 7: 10. Plot of the values of rj that reproduce the observed rjB, as a function of mi 
for p = 1. The different curves correspond to = 0.5 (dotted), yjs = 1 (dashed), yp = 1.3 
(long-dashed), yp = 1.7 (dot-dashed) and y^ = 2 (solid). The gray band is the indicative 
optimal range of values of = 5 x 10"'^ 5 x 10~^. 

values of < 5 X 10^^ ^ 5 x 10^^ are needed. Therefore, in these examples, the observed 
asymmetry is reproduced for reasonable values of ~ 0.01 only for yp = 1. For this 
value there is indeed a compensation between very large values of the CP asymmetries, 
\ei\ ~ 10^^ for 1] ~ 10~^, and an additional wash-out suppression ~ 10~^ coming from 
AL = 2 processes so that the resulting efficiency factor ~ 10^^. 

We have again summarized the situation plotting in Fig. [3, the value of rj such that 
\Vb\ = rj^'^^ as a function of mi. We show five examples for yis = 0.5, 1, 1.3, 1.7, 2. One 
can see that this time 1] naturally falls in the optimal range rj = 5 x 10^^ ^ 5 x 10~^ (the 
grey band) only for for = 1^1.7, therefore requiring some amount of tuning. Moreover 
this occurs not for all values of mi and in particular values y/3 ~ 1 require m; > 0.1 eV. 
Therefore, an improvement of the upper bound on mi (cf. (1571) ) will enforce higher values 
of yi3 corresponding to higher values of the RH neutrino masses that in turn imply higher 
values of the initial temperature of the radiation dominated regime that, in an infiationary 
language, corresponds to the reheating temperature of the Universe. 

4.4 Reheating temperature constraints 

The model predicts a RH neutrino spectrum with a mild hierarchy and an overall scale of 
values of the RH neutrino masses that is quite large. For NO one has M3 ~ 10^^ GeVy^, 
where M3 is the mass of the lightest RH neutrino. For 10 the scale can be about five 
times lower if one considers quasi-degenerate light neutrinos with masses mt > 0.1 eV, 
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close to the current cosmological upper bound (cf. (157|) ). 

Since RH neutrinos are produced by thermal processes, this in turn implies a lower 
bound on the reheating temperature given approximately by [5] Treh ^ ^3/[^_b(-^3) — 2] ~ 
10^^7/|GeV in the case of NO and about five times lower for 10. Are such high values 
of the Treh possible ? Since our model is supersymmetric, the well known upper bound 
^reh ^ lO^^^'^GeV from the avoidance of the gravitino problem potentially applies |48j . 
It is clear that such low values cannot be obtained in the presented version of the model 
since, even taking the low value yp ~ 0.1, one obtains Treh ^ 10^^ GeV for NO. 

There are two possible kinds of ways out remaining within thermal leptogenesis. The 
first one would be to circumvent the gravitino problem and indeed a few solutions have 
been proposed ^9]. The second kind would be to modify the model in such a way that the 
overall RH neutrino mass scale is lowered. It would not be difficult to envisage different 
schemes to this extent. 

On the other hand, in the case of 10 the second strategy cannot be invoked for the 
simple reason that large values of the RH neutrino masses Mi > 10^^ GeV are neces- 
sary in order to get a strong additional wash-out suppression from AT = 2 processes to 
compensate the very large values of the CP asymmetries. In other words large reheating 
temperatures Treh ^ 5 x 10^^ GeV cannot be avoided in the case of 10. 

5 Conclusions 

A typical outcome of extensions of the standard model that attempt to explain the features 
of the lepton mass spectrum on the basis of fiavor symmetries is that small quantities such 
as the charged lepton mass ratios, Am^^^/Am^j^, 6'23 — 7r/4 and ^^13 are proportional to 
some power of small symmetry breaking parameters rj. By keeping only the leading order 
power of 7], in some cases the number of independent parameters becomes small and 
these models can be rather predictive, with characteristic relations among the observable 
quantities. Well known examples are relations between the neutrino oscillation parameters 
and the branching ratios of lepton fiavor violating processes, such as /i — 67, r — /i7 
and T 67. It would be quite interesting to include in this kinships also the baryon 
asymmetry, which, in the context of leptogenesis, is naturally related to lepton masses 
and mixing angles. 

With such scope in our mind, we have discussed the constraints on the leptogenesis CP 
asymmetries in models possessing a fiavor symmetry. We have derived general conditions 
for the vanishing of the CP asymmetries in the limit of exact fiavor symmetry. We 
have shown that, if the three RH neutrinos belong to an irreducible representation of 
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the flavor symmetry group, then the total CP asymmetries are zero in the hmit of exact 
flavor symmetry. More precisely, for non-degenerate RH neutrino masses, the total CP 
asymmetries are of order r)^ and the flavored ones ej-^a are of order r). If the RH neutrinos 
are not in an irreducible representation of the flavor group, we have derived a necessary 
and sufficient condition for the vanishing of the total CP asymmetry in the symmetric 
limit and we have discussed it in several particular cases. For instance, if the action of 
the symmetry on RH neutrinos is Abelian, then in most cases the CP asymmetries are of 
order 77° and we should invoke additional washout suppression to reproduce the observed 
baryon asymmetry. 

One interesting example of vanishing leading-order CP asymmetries is that of a model 
symmetric under x Z3 x U{1)fn, built to reproduce tri-bimaximal Icpton mixing. In 
this model RH neutrinos arc in a triplet of A4 and = 0{7f). The model is rather 
constrained. Once the parameters arc fixed to match Aml^i and Am^^^, there is only 
one relevant phase (/?, which can be thought of as a function of the lightest neutrino mass 
mi. The RH neutrino spectrum depends only on an additional, 0(1), parameter yg. Both 
normal and inverted neutrino mass ordering can be reproduced. For normal ordering 
nil is essentially fixed in a small range around 0.005 eV and the phase should be very 
small. For inverted hierarchy there is a lower bound on mi of approximately 0.017 eV 
and there is much more freedom for the phase. Given this rather constrained framework 
it is not guaranteed that a successful leptogenesis can take place at all and we computed 
the washout effects and the resulting baryon asymmetry. The dynamics of the model 
is quite interesting since the RH neutrinos have similar masses and they all participate 
to generate the baryon asymmetry. At the same time, to a good approximation, the 
asymmetry produced in the decay of one RH neutrino is not washed-out by the other 
heavy neutrino inverse processes, since the interactions of the three RH neutrinos with 
the light leptons are almost orthogonal to each other. For normal hierarchy we find 
that the observed baryon asymmetry is reproduced for values of the symmetry breaking 
parameter 77 in the range 0.005 -j- 0.05, which nearly coincides with the natural expected 
range in this model. This prediction is rather stable with respect to variation of y^j. For 
inverted hierarchy we find solutions in the parameter space, but they are less stable. If 
0.005 < 7] < 0.05, there is only a rather limited range of allowed values for y^. On the one 
hand, for small values of the baryon asymmetry is typically enhanced, compared to 
the normal ordering case and we cannot go below ~ 1, with a reheating temperature 
not lower than 5 x 10^^ GeV. On the other hand, as soon as y^j exceeds 2 the suppression 
from the washout becomes huge and the baryon asymmetry goes rapidly to zero. 

It is interesting that though the 10 case has a much wider parameter freedom compared 
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to the NO case, it certainly appears less attractive from a cosmological point of view. 
Future improvements on the measurements of the absolute neutrino mass scale will allow 
to test our results. 
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